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Forced advection of passive tracer, 8, in nonlinear relaxational medium by large scale (Batchelor 
problem) incompressible velocity field at scales less than the correlation length of the flow and larger 
than the diffusion scale is considered. Effective theory explaining small scale scalar fluctuations is 
proven to be linear, asymptotic free (downscales from the scale of the pumping) and universal. Only 
three parameters are required to decribe exhaustively the small scale statistics of scalar difference: 
two velocity-dependent ones, average and dispersion (A and A respectively) of the exponential 
stretching rate of a trial line element, and a, standing for average rate of linear damping of small 
00 ' scale scalar fluctuations, a is an explicit functional of potential chracterized medium nonlinearity 

0^ ' and amplitude of 8 2 flux pumped into the system. Structure functions show an extremely anomalous, 

intermittent behavior: (|<5# r | 9 ) ~ r ^ q , C? — m i n ^<Zi \J [s] + — SJ ' ^° dissipative anomaly is 
f> | found in the problem. 

m _ 

I. INTRODUCTION 



Turbulence is very nonequilibrium state of nature, which becomes stationary if energy is supplied permanently at 
large scales. To construct a theory of turbulence means to describe temporal and spatial distributions of velocity 
and of variety of different thermodynamic characteristics of the fluid, i.e. density, if turbulence is compressible, 
temperature, if thermo-advection is applied, relative concentration of components in case of multi-component (color) 
flow, magnetic field distribution in a conducting fluid etc. Dynamics of different fields describing a real turbulent flow 
is both nonlocal and nonlinear. We call the general situation active to emphasize the reciprocal character of interaction 
between velocity field and thermodynamic characteristic(s). However, sometimes the effect of a thermodynamic field 
on the velocity distribution is suppressed. It takes place , for example, if scales are separated: a typical spatia-temporal 
scale of velocity is much larger than one of a thermodynamic quantity. The case, when it is theoretically justified 
to neglect the effect of back reaction of thermodynamic field on velocity field in comparison with ones of advection 
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and nonlinearity is called passive. The passiveness does not necessarily means linearity. Moreover, our objective is to 
' , study passive yet nonlinear situation. 

We consider dynamics of a thermodynamic quantity 9, 

H{9} = J dr[^(V8) 2 + U(9)], (2) 

b : 

where H{9} is a positive definite thermodynamic functional of the system, U{6) is confined ( U — > +oo at 9 — > ±oo) 
potential, k is diffusion coefficient ; (f>(t; r) stands for statistically steady forcing providing a constant supply of 
otherwise relaxational ^-dynamics at large scales []. 9 is imbedded in a turbulent flow, i.e. the temporal derivative is 
extended by sweeping term 

d d 

where incompressible velocity field u(i; r) is prescribed to be known statistically. 

We aim at finding the statistics of the passive scalar 9 fixed by (|f[]|) in the inertial interval of scales i.e. for scales 
that are less than both the velocity correlation scale, L u , and the scale of the scalar supply, L, and larger than the 



*based on a talk delivered at "International Conference on Turbulence: Challenges for the 2fst Century" held at Los Alamos, 
May 18-21, 1998 and dedicated to R.H. Kraichnan 70th birthday. 

x We will discuss here the simplest case possible, when the thermodynamic field is a scalar. Notice, however, that generalization 
of the discussed theory for a vector or generally tensorial object is possible. 



1 



diffusion scale. Incompressible velocity field at those scales is modeled by the first term of its local expansion in the 
radius vector connecting a reference point with the current one, 



u(t;r)=*(t)r, (4) 

where a(t) is d x d traceless random matrix of velocity's derivatives. 

The problem (Q-Q) describes forced advection of a scalar pollutant in the viscous-convective range absorbed or 
generated depending on the sign of the nonlinear rate dgU(9), for example via a chemical reaction with other species 
presented in abundance in the flow. The problem is of fundamental importance for geophysical atmospheric turbulence 
(see m for review) . Other relevant phenomenon is turbulent thermo-advection in a cell attached to thermal bath (see 
H and reference therein). Then, dgU(9) is nonlinear heat transfer coefficient and 9(t; r) measures local deviation from 
the bath temperature. Many regimes of premixed turbulent combustion are also governed by (|l^§|) @,. The last (but 
not the least) situation to mention is a phase ordering in a system described by scalar nonconserved order parameter 
(a very well known object of the phase transition theory, see [fij-fl for reviews) advected by large scale turbulent flow 
(see also || on some discussion of reciprocal description of advection and critical dynamics). 

Our consideration will be based essentially on understanding, results and general terminology emerged from studies 
of the pure problem of passive scalar advection (no medium effect at all, U = 0) having an almost five decades of 
history (see Obukhov and Corrsin papers ||, for the earliest contributions). Batchelor |TlJ has pioneered study 
of the smooth velocity field limit (Q), which nowadays has grown to be (through important contributions of many 
people fiHH) one of the most advanced theory in the field . Temporal short-correlated but spatially non-smooth 



model of velocity, one which gave more than two decades later the first ever analytical evidence of intermittency in 
turbulence was invented by Kraichnan. Structure functions of scalar difference in the convective range, 

S q (r)={\6(t;r)-6(t;0)\' l )~r 6 '>, (5) 

became the key object in the intermittency study. The anomalous scaling, A2„ = — £,2m describing the law of the 
algebraic growth with L/r of the dimensionless ratio, S2 n (r)/ [^(r)]™, was shown to exist generically p4j--p6|. The 
anomalous exponents were calculated perturbatively in expansions about three non- anomalous ( A2 n = 0) limits, of 



large space dimensionality d [ p4|j27| , of extremely non-smooth [E5 28 1 and almost smooth |E6[ velocities respectively. 
A strong anomalous scaling (saturation of ^2n to a constant) was found for the Kraichnan model at the largest n 
by a steepest descent formalism p9| , ^0[ . Although the restricted asymptotic information about anomalous exponents 
in the model is available a future possibility to establish rigorously complete dependence of £2™ on n, d and degree 
of velocity non-smoothness seems very unlikely (in a sense, recent Lagrangian numerics |]3l|| compensates the lack of 
rigorous information). 

The problem (0-^) discussed in the present paper is also showing anomalous scaling, £2™ < n^2, which is resolved 
analytically for any kind of temporal correlations in a and an arbitrary (yet confined at 9 — » ±00) potential U{9). 
An important circumstance helping to find the general answer is scale separation between scalar fluctuations at the 
integral scale, L, and the current one, r, from convective range. For general U(9), S9 r obeys the same statistics as one 
would expect from an auxiliary (linear!) problem with quadratic potential, U*(9) — a9 2 , where a is given by average 
of dgU (9) with respect to single point scalar distribution, V\ ~ exp [— U(6)/xo], with xo = dt (</>(t; 0)0(0; 0)). a 
is always positive, i.e. at the smallest scales the effect of nonlinearity, generally alternating between damping and 
acceleration, is reduced to a pure linear damping. Linear problem, by itself, happens to be solvable not only in the 



case of the short-correlated velocity 32 but for a(t) statistics of a general position (providing correlation functions of 
a(t) taken at different moments of time decays with the time shift faster than algebraically). Finally, we have found 
anomalous exponents 



£q = min { q, 



2aq A 
~~A A 



(0) 



where A and A are respectively average and dispersion (with respect to a averaging) of the rate of line stretching, 
Aft) = t^ 1 \n[R(t) / R(0)] , with R(t), satisfied to R(i) = &(t)R(t). To make the statement we used the method of 



The anomalous behavior in the model under consideration differs from one perceived in the Kraichnan model. First 
of all, ^2n as a function of n does not saturate to a constant at the largest n but keeps growing with n as yfn. Second 
( and major) difference is associated with the concept of dissipative anomaly. It is generally accepted to talk about 
dissipative anomaly if some object calculated at zero dissipation ( k = 0) does not coincide with its K — > counterpart. 
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In the Kraichnan model the anomalous scaling coexists with dissipative anomaly |23j]. However, the nonlinear problem 
(|l[^), as well as its linear descendant, shows no dissipative anomaly while the anomalous scaling is present. We base 
the important conclusion on the following no anomaly criterium (which, we believe, is general): if zero dissipation 
analysis produces normalizable and everywhere positive solution for the Probability Density Functional (PDF) of 
fluctuated field ( 9 in our case) then the dissipative anomaly is absent. 

The problem is formulated in Section II. To describe the scalar fluctuations at a current scale from convective range 
we show how to integrate out the large scale contribution in Section III. The scale separation results in suppression of 
nonlinearity. The effective small scale theory appears to be a linear one with uniform damping. All the final answers 
emerging from study of the linear problem are presented in Section IV. The last Section V is reserved for conclusions. 



II. FORMULATION OF THE PROBLEM 

(|l|||) describe advection of a passive scalar 9(t;r) by the smooth incompressible velocity field (0). The scalar is 
forced by random field cj)(t; r), which for a sake of simplicity is considered to be Gaussian thus fixed unambiguously 
by 

(0(ti; n)0 (t 2 ; r 2 )) = X (|n - r 2 |)<5 (t a - t 2 ) , (7) 

where the function X (r) decays fast enough if r exceeds the integral scale L < L u . xo = x(0) is the fhrx of 9 2 pumped 
into the system, a is a random in time matrix process described by its PDF, <f>{ a (t)}, which is supposed to be known. 
Diffusion is supposed to be small, such that the range of scales in between = v / k/[S'/t] 1 / 4 ( S and r are typical 
values of the strain and velocity correlaion time respectively) and L is sufficiently large, L/rd ^> 1. 
We will be mainly aiming to find the two-point scalar PDF, 

V 2 (x+, z_|r) = (<S (z_ - 9{t; 0) + 9{t; r)) S (x+ - 9{t; 0) - 9{t; r))) , (8) 

and the scalar structure functions, 

S 2n (r)^([9(tr)-9(t;0)] 2n ), (9) 

where averaging, (•••), with respect to both a{t) and cf>(t; r), is assumed. 

Other important objects used in the course of the forthcoming calculations will be the two point scalar PDF, 
measured for particular o(t) configuration , 

Qi (xi, x 2 \r h2 ; t; {a{t'); -co < t' < t}) = (S (a* - 9{t; n)) 6 (x 2 - 9(t; r 2 ))) , (10) 

and the single point scalar PDF 

V 1 (x) = {6(x_-e(t;v))). (11) 

Notice, that the last object does not depend on the velocity field statistics because of spatial homogeneity assumed. 
Deep inside the convective range (at L ^> r\ 2 ), \9\ — 9 2 \ <C \9\ + 9 2 \, and ( px| ) can be decomposed into the product 

G2 (9i,9 2 |n, 2 ; t; {&(t');-oo <t'<t}) = Vx(6x) * Q- {9 X - 9 2 |n - r 2 ; t; {&(t');-oo <t'<t}). (12) 

The average of (|l^) over a reads as 

Vi{e u e 3 \ri2)=Vi(9 l )*V-(0i-B2\rx2). (13) 

The assumption on the absence of the dissipative anomaly in the case of a very small diffusion lies in the core of 
our consideration. The formal consequence of the statement is the possibility to drop off the dissipative K-dependent 
term from ^ already on the dynamical (yet unaveraged) level. The no-anomaly assumption will be justified by the 
positivity and normalizability of the derived answers for PDFs. 
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III. REDUCTION OF THE NONLINEAR PROBLEM TO A LINEAR ONE 



In the absence of diffusion can be integrated along the Lagrangian trajectories (characteristics) 



dt: 



7 0(t';p(t')) 



du 
~~dJ 



e(t'-p(t')) +<p(t';p(t')) 



—p(t')=a(t')p(t>), p{t) 



oo < t' < t. 



(14) 
(15) 



Notice, that the nonlinearity leads to scalar generation in the region of convex (dglf > 0) potential while it dumps 
scalar fluctuations if dell < 0. Fokker-Planck equations (see Q for similar calculations) derived out of ( |l4| , |l5| ) by 
means of direct averaging over the Gaussian noise, <f>, are 



_ dU (9) „ 2 



dU ((h 



i,i=l,2 



Vi = 0, 

^2=0, 



(16) 
(17) 



where G2 is not stationary, since it does dependent on time explicitly through a(t). Integrating ( p"7| ) with respect to 
9+ = 61 + 02 and assuming that the integral is formed at \0i — Qi\ <^ \6\ + 62], where ([[2]) is valid, we arrive at the 
close equation for the scalar difference PDF, 



{d t + (a^(t)r^ - ad x x) - 2 [ X (0) - X (r)] d 2 x ) Q- (x \r;t; {a(t'); -00 < t' < *}) = 0, 
where a is defined as the following average over the large scale 6 statistics 

oo 

'd 2 U{9)\ f ja d 2 U{0). 



d9 2 



LS 



d6- 



de 2 



-Vi (9) . 



(18) 



(19) 



The normalized and everywhere positive solution of (|T^) is 

exp[-U(6)/ X o] 



Pt (9) 



J d9eiq>[-U(6)/xo} 



(20) 



Substitution of (M) into (h9J) gives 



dU(9) 
d6 



I d9 



LS 



dU(9) 



exp[-U(0)/ X o] 



Xo 



Xo f d9cxp[-U(9)/xo} 



(21) 



i.e. a is principally positive constant, does not matter what is the form of the potential U{9) be (provided it grows 
with \9\ at 9 — > ±00). Therefore, we have found that at the smallest scales regions of scalar generation are suppressed 
statistically. 

On the basis of ( |l6| ) and ([l8]) we conclude that, from the point of view of the small scale statistics of scalar difference 
our problem is equivalent to the linear one, with dU(9)/d9 being replaced just by a9. In other terms, we may proceed 
averaging the linear dynamical equation 



d t 9 + a" u {t)r^9 = -a>9 + <j>{t; r), 



(22) 



instead of the original nonlinear one. The steady distribution of the scalar difference underlined (|22j) was the subject of 
the recent paper [[32| , the method and results of which will be briefed and generalized (for the case of finite correlated 
velocity) in the next Section. 
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IV. VELOCITY AVERAGING. ANOMALOUS SCALING. 



The linear analog of ( [hi] ) is 



0{t; r) = / d£ exp [-at'} <t> (*'; p(i - t')) 
o 



(23) 



For the purpose of the 2n-th structure function calculation it is utmost enough to consider the simultaneous product, 
F\...2n = - ■ ■ #2n)> which is according to (|23|]15|) is 



{l,-,2n} I n °° 

e (ny^ e_ 

{*!,— i»2n} \fc=l 



= Texp 



dt'a(i') 



dW(i) 



di 



= *(f)W(f). 



(24) 



(25) 



Calculation of F\...2 n is essentially simplified for the collinear configuration, r, = n fj, when the 2n x (d— 1) parametric 
average (p4h is reduced to the following single-parametric one 



{l,-,2n} / n 

e (n 

{*!,... ,»an} \*=1 o 



dt k e~ atk 



(26) 



with Tij = |rj — r,-|. Here, the longitudinal stretching rate, r/(t) = hx\W(t)n\, is the only fluctuating quantity left. 
The a = version of was calculated in jlT]] for the d = 2 case and generalized for any d > 2 in [jl8| via a change 
of variables and further straightforward transformation of the path integral standing for the average over a(t). It is 
shown that at the largest times ^-measure is a shifted Gaussian one, 



T>rj{t) exp 



dt 



2A 



(27) 



characterized by two parameters only: average, A, and dispersion, A, of the Lyapunov exponent, J dt'r](t')/t. ( p7| ) 

r— ■ 

applied to (£6|) produces 



F\...2n 



III 



JJ dtidrji 



exp 



A 



A 



2 -, {1,— ,2n} „ 
{As 4 ,- 1 fc a „}t=l 



where r/i ( i < n) integrations are not restricted, < £ n < • • • < ti < oo, t n+ i = rj n+ i = 0, ti± 
equivalent notations for rj) and 



exp 



2A/ 



V2?rAi 



(28) 
— £fe (with 

(29) 



Both A and A are unambiguously fixed by <f>{a(t)}. The integrand of (pq ) decays exponentially in time with the 
major contribution into the integral formed at ti ~ 1/a. The leading term does not depend on any ry and gives no 
contribution into 2n-th order structure function. The first actual r-dependent contribution stems from n—1 temporal 
integrals formed at t ~ and one at U ~ t t ~ In [L/r] / max{a, £>}. This special integration brings a spatial 
dependence into the object, therefore on a single distance. Generally, there exists a variety of terms with all the 
possible combinations, like term with k integration formed at r, while n — k ones at r r , and therefore dependent 
explicitly on 2(n — k) spatial points. However, we are looking exclusively for a term dependent on all the 2n points 
since only such a term contributes S^nM- ^ is really simple to calculate the scaling of this term making use of the 
temporal separation, r r ^> r. Indeed, the large time contribution may be extracted out of pi) in a saddle-point 
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calculation. Variation of all the exponential terms in (|2§|) with respect to U gives a chain of saddle equations. The % 
functions in the integrand of (Eq) limits the ij integrations from above by In [L/r] . Therefore, the desirable 2n-points 



contribution forms at ti = y A/ [A (2cmA + A 2 /2)] In [L/r], and rji — In [L/r], where it is assumed that in the leading 
logarithmic order there is no need to distinguish between contributions of different separations . Substituting the 
saddle-point values of ti and rji into ( |28| ) one arrives at the anomalous part of (^|), with q = In. The normal-scaling 
counterpart of (^) originates from expansion of x( r ) from (^8|) in a regular series in r 2 . 

The basic physics of nonzero £2™ (means deviating from the naive balance of pumping and advection) and generally 
anomalous (£2™ < ^£2) scaling at a > can be stated quite clearly. According to (|2^) the advection changes scales 
but not amplitude, while the amplitude of injected scalar field decays exponentially from the time of injection at the 
constant rate a. The temporal integrals in (^8|) forms at the mean time to reach a scale which is proportional to the 
negative log of the scale. However, the effective spread in the factor by which amplitude has decayed, upon reaching 
a given scale, increases as scale decreases. It is why £2™ > 0. Also there is more room for fluctuations about the mean 
time due to the interference between the exponential decay of the scalar amplitude and fluctuations of the stretching 
rate rj. Thus intermittency increases with a scale size decrease. 

Another way to derive ( ||) out of (^3|) is to construct S2n(r) directly. It is easy to check that the structure functions 
of different orders are produced by the PDF satisfied to 

Xr^'^dy+^l^drV^ + ad x (xV-) + [ X o - X (r)] r 2 d 2 x V- = 0, (30) 

The solution of (|3^), in the regime where you can neglect the x _ dependent term is 

+ +ioo 



P_ (jB |r) = JL-L / „„ 



s + 1 



r 



L 



y/d 2 /4+as/[D(d~l)]-d/2 



a s . (31) 



0+— ioo 



Here, a s is a function fixed by matching at the integral scale, roughly, V- (x \L) ~ V\ (x), where V\ (x) is given by 
(pp|). The PDF ( |3l| ) is positive and normalizable, that, therefore, confirms the initial hypothesis on the absence of 
dissipative anomaly. Also, (|3l| ) shows that (||) holds for general (not only even integer) positive q. 



V. CONCLUSION 

We have shown that the nonlinear problem (fi]^|) is reduced to a linear one at the smallest (still from convective, not 
dissipative, range) scales. The asymptotic theory remembers about initial nonlinearity through the effective damping 
coefficient (p^). The linear problem was solved for the general case of arbitrary correlated in time large scale velocity 
field. 

The most important feature of the problem appears to be the absence of dissipative anomaly, the point which was 
guessed initially. Selfconsistency of the hypothesis was confirmed afterwards by checking positivity and normalizability 
of the final expression (|3l]) for PDF. Of course, the absence of dissipative anomaly is, by no means, a common situation 
in turbulence. It is however suggestive to start analyzing any new turbulent problem from the simple "no anomaly" 
test. 

We discussed only the Batchelor case of large scale velocity. The simplification appears to be a very important both 
for the fact of absence of dissipative anomaly and solvability of the problem. The Batchelor case is very special, since 
the Lagrangian dynamics of n particles, generally described by n(d — I) degrees of freedom, is reduced to dynamics 
of d — 1 eigenvalues of stretching matrix. This lies in the core of the Batchelor problem's solvability. Also, in the 
Batchelor case scaling dimension of eddy diffusivity operator coincides with one of the a (damping) dependent term. 
The coincidence of exponents explains the anomalous scaling, particularly the continuous dependance of the exponents 
on a. Any multiscale velocity field (say taken from the Kraichnan model) leads, first, to appearance of the dissipative 
anomaly already on the medium free ( U = 0) level, and second to dis-balance of the scaling dimensions of advective 
and medium-originated contribution into the eddy-diffusivity operator, resulting in complete screening of any medium 
effect in the convective range. We conclude by this guess, which is rather brave (and , of course, is not rigorous at 
all). More studies, first of all on the nature of dissipative anomaly, are required in this direction. 

I thank L. Kadanoff, I. Kolokolov, B. Meerson, A. Patashinski, R. Pierrehumbert, A. Polyakov, B. Shraiman and 
V. Yakhot for inspiring discussions. Valuable comments of P. Constantin, G. Falkovich, R. Kraichnan, V. Lebedev, 
and M. Vergassola are greatly appreciated. The work was supported by a R.H. Dicke fellowship. 
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